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Abstract 



We derive a formula for computing the transition rate for a process involv- 
ing particles with momentum much higher than the temperature and chemical 
potentials in a plasma by using an effective field theory approach. We apply 
it to collision of charged particles with hard momentum inside a QED plasma. 
The Debye screening effect and the damping of a charged particle moving in 
QED plasma are studied. Using the Bloch-Nordsieck resummation, the in- 
frared divergences due to the absence of magnetic screening for QED plasma 
are shown not to appear in physically measurable rates. The soft plasmon 
absorption and emission for charged particles are discussed. 
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I. INTRODUCTION 



Effective field tlieory lias been proved to be powerful in dealing with different energy 
scales. Consider a process that contains particles with hard momentum (hard momentum 
means one of the component of the four momentum is larger than the temperature and the 
chemical potentials.) happening inside of a plasma. There are two relevant energy scales: 
the hard momentum scale and the soft momentum scale which is the thermal scale, i. e. order 
of the temperature T or the chemical potential. It is convenient to first "integrate out" the 
hard field degree of freedom with the momentum scale much higher than the temperature 
and the chemical potentials and end up with an effective action for the soft field degree of 
freedom. Taking the thermal average of the soft field degree of freedom yields the transition 
rate for the process. In this way, we can deal with the hard momentum scale and the thermal 
scale separately. Applying this formula to the hard process involving charged particle inside 
a QED plasma automatically incorporates the Debye screening effect and enables us to study 
the infrared physics for QED plasma due to the exchange of soft photon between the charged 
particles and the plasma, such as the damping of a charged particle moving in QED plasma 
as well as the emission and absorption of soft plasmon. 

It is well known that the infrared divergences coming from the radiative correction due 
to the virtue soft photon are canceled by the real soft photon emission |I|J^. When the 
system is heated to a finite temperature, it is even more sensitive to the infrared region for 
the photons because of the bosonic distribution factor n{uj) ~ T/u. It is interesting then to 
examine the cancellation of the infrared divergences for QED at finite temperature where 
a large population of thermal soft photons due to the Bose distribution enhances both the 
emission and the absorption of the soft photons for the charged particles. In the literature, 
there has been study of this problem for the free thermal photon background 0,^. Our 
formula enables us to provide an analysis of these emission and absorption processes for the 
interacting photons inside QED plasma. The interaction between photons and the plasma 
brings the Debye screening mechanism to shield the electric force but not the magnetic force. 
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Therefore, potential infrared problem could occur due to the exchange of soft magnetic 
photons. In fact, there have been many discussions on the infrared divergence for the 
damping rate of a charged particle moving in QED plasma P . As discussed later in section 
p.11 C| , this arises because the unscreened magnetic force does not allow the charged particle 
moving with a definite momentum. For a physically measurable damping rate, this infrared 
divergence is cut off by the momentum resolution of our experiment. In general, we show 
that the infrared divergences due to the unscreened magnetic force cancels out by using 
the Bloch-Nordsieck resummation. We also discuss the correction to the transition rate due 
to the soft plasmon emission and absorption. At the massless limit of the plasmon, we 
reproduce the previous result 0. 

The organization of this paper is as following. In next section, we develop first a general 
formula for calculating the transition rate for a hard process in a thermal background. In 
section III, our formula is applied to process inside QED plasma to study the Debye screening 
effect, the damping of a charged particle moving inside QED plasma, and infrared physics 
of QED plasma. Conclusions are drawn in section IV. 

II. THERMAL AVERAGE OF TRANSITION RATE 

In this section, we derive a general formula for computing the transition rate for a process 
involving hard momentum particles with the appearance of a relatively low (compared with 
the hard momentum scale) temperature T thermal background. The interaction between 
the thermal background and the hard particles will be considered. 

Let us denote the fields with hard momentum by 0h and the fields with soft momentum 
by 0s- Since, we consider the case where the temperature is much lower than the hard 
momentum, for the hard momentum fields, it is valid to treat them as at zero temperature. 
The thermal averaged transition rates reads 

(R) = ;^^EK/,^K,^)l'e-^^''""^'"''^'^"^^ (2-1) 

t m,n 
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where Tt is the total time for the interaction to take place and i, f label the initial and final 
states of the hard momentum particles respectively, m, n represent the energy eigenstates 
generated by the soft momentum fields (ps- A thermal distribution factor has been included 
with P = 1/T, fii are the chemical potentials corresponding to conserved charges Qi. En 
and Qi{n) are the eigenvalues of the Hamiltonian and the charge Qi for state \n). Z is the 
partition function for the soft momentum fields 

Z = Tr (e~^^^~^^ ^''5')) = J2 e"'^^^""^" /^.Q.W) _ (2.2) 

n 

To be more specific, we shall consider initial and nf final hard momentum particles 
are involved. Let us use p'j- with / = l,2,---,n/ and pi with i = 1, 2, ■ ■ ■ , denote the 
momentum of the final and initial particles. Therefore, 

(2.3) 

where flf are the phase space volume within which the detectors are hunting for the final 
products. V is the total volume of the spatial box in which the wave function of the initial 
and final states are normalized. Expression ( p.3|) has following structure: matrix element 
{m,p[,p'2, ■ ■ ■ ,p'n j:\pi, P2, ■ ■ ■ ,Pni,n) iuvolves time evolution from far past to far future while 
its complex conjugate evolves from far future to far past; the thermal distribution factor 
exp{— /?(£'„ — J2i f^iQi{n))} causes an evolution in imaginary time by amount ip. A real 
time formalism for the thermal field theory can conveniently incorporate this structure 
so that we can express it in terms of functional integrals defined along paths in complex 
time plane. Slicing the time evolutions into pieces by inserting complete sets of states as we 
usually do to derive the path integral formalism and using the reduction formula to write 
the matrix element in terms of Green's function, Eq. ( |2.3| ) may be expressed as 



X / ^^0h n / d^^f9*f{p'f,Xf)-dxo4>h{xf) n / d^^iMPi,Xi)-dxo4>h{xi 

JPl f=l ^ i=l ^ * 



X exp{iSp[(l)s] +2S'pj0h,0s] +i5'p2[0h,0s]} • (2.4) 

Here the notations need to be explained. fi{pi, Xi) and gf{pf, xj) are the wave functions for 
the initial and final single particle states respectively. The path integrals are along contours 
in complex time denoted as r plane: contour Pi starts from r=— oo to r=+oo— ie; contour P2 
starts from r=+oo— ze to T=—oo—2ie; contour P includes Pi, P2, and another piece from 
r=— 00— 2ze to t=— 00— i(3. Boundary conditions are 0s(t=— oo,x)= ± 0s(r=— 00— 2/5, x) 
for the soft bosonic (+) and fermionic (— ) momentum fields while 0h(T= ± 00, x) for 
hard momentum fields. Actions Sp[(j)s], SpJ^h, 0s], and S'pJ^h, 0s] are defined with the time 
integral being the contour integral in the complex time r plane along P, Pi, P2 respectively. 
S'p[0s] represents the part of the action contains only 0s while 5'pj0h, 0s] and 5'p2[0h,0s] 
contains both the part containing only 0h and the interaction terms between 0h and 0s. The 
time component of vectors Xi,Xf belong to path Pi while the time component of vectors x'^, x'^ 
belong to path P2. This type of path integral is used usually in the real time formalism for 
thermal field theory 

From an effective field theory point of view, we can first integrate out the hard momentum 
degree of freedom. The formula above may be rewritten as 

x{Pl,P2,■■■,Pn^p'l,P2,■■■,P'nf)'f>s,P2 (2.5) 

or more compactly 

{R) = ^ [ I^0sexp{z5p[0s]} {m^MMf)<t>.,P. , (2.6) 
if jp 

where {f\i)(f,s,Pi represents the matrix element between the hard momentum states \i) and 
(/I with the background soft momentum field 0s appearing. Subscripts Pi,P2 denote the 
time argument in field 0s belonging to contours Pi, P2 respectively. We note that the thermal 
effect has been taken into account by the "strange" path integral that we use and the chemical 
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potential for the conserved charges are incorporated by including them in the action 5'p[0s] 
as the Lagrangian multipliers. It is convenient to introduce source terms for field with 
time argument belonging to paths Pi and P2 respectively and define following generating 
functional 

^'"^^'^''^^ = I exp ^tSp[<P,]+t J^dr J d^Ur, x)03(r, x)+2 J^dr'J d^'j^ir', x')0:(r', x') 

(2.7) 



With this, 



{R) = ^ Af, 



l_5_ 



A* 



if 



1_5_ 

i 5]l 



where the functional ^/j[0s] is defined as 



(2.8) 



(2.9) 



When the argument for Afi is a functional derivative operator, the power series expansion 
of the functional is understood. We note that the path integral along contour P is simply 
another way to write the thermal average. Generating functional ( ^.7[ ) may be written as a 
thermal average in an operator form: 



-Tr 



f3{H- t^rQr)rp Li J^J^ J ^ixi, (t ,x)^, (i ,x) \ f-i j'^Jt'^ dx'je (i' ,x' )<^| (x' ,i' ) 



(2.10) 



where T4. means time-ordered product according to t and T_ means anti-time-ordered prod- 
uct according to t' . Here the Hamiltonian H and the charges Qi are operators for the degree 
of freedom corresponding to field operator 0s- 

The formula above enables us to deal with the two scales involved in the process sepa- 
rately. The characteristic energy scale for amplitude {f\i),p^^Pi is the hard momentum scale. 
We can then treat (p^^ as a slowly variant background field which enables us to use, for ex- 
ample, the semi-classical method. After we completed the hard momentum scale, we only 
need to deal with the soft momentum scale physics left. 



We note that it does not require the momentum transfer to be big to get the formula 
above. All we need is that the particles involved has high momentum so that the thermal 
bath does not contain a noticeable amount of these hard momentum particles which we 
want to scatter and detect. This formula is also valid for the case where we just formed 
a plasma which is weakly coupled to another kind of particle, for example, the neutrinos. 
The plasma does not have enough time to produce the thermo-neutrinos, so that we can 
treat the neutrino fields the same way as 0h without requiring neutrinos to carry only hard 
momentum because the neutrino field is effectively at zero temperature. Practically, we can 
do a perturbation in the weak coupling constant to evaluate {f\i),j>^ for using this formulation. 

A concrete formulation for computing non-relativistic reaction rate has been well devel- 
oped in reference 

III. PROCESS FOR HARD CHARGED PARTICLES IN QED PLASMA 

We now apply the formulation in previous section to study explicitly a collision involving 
hard charged particles with the appearance of QED plasma. Let the process start with 
initial particles and end up with Uf final particles plus soft photons. We assume that the 
initial particles have the four-momentum pi, with i = 1,2, ■ ■ ■,nj and we detect the final 
particles with momentum p'j with / = 1, 2, ■ ■ ■ , rij but left the soft photons undetected. 
In view of the formula in the previous section, 0s now represents the thermal photon and 
thermal electron fields while 0h represents the hard charged particles. Integrating out the 
degree of freedom corresponding to 0h to get an effective functional integrand for 0s means 
to get the zero temperature transition rate for hard charged particles with the appearance 
of the thermal QED background field. Expression ( p.6|) for transition rate reads 
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i=l f=l 

xG[xi, X2, ■ ■ ■ , Xn,, 2/1, 2/2, ■ ■ ■ , Vrif] Ap^]G* [x[, X2, ■ ■ ■ , x'^^, y[, y'^, ■ ■ ■ , y'^y, Ap^] . (3.1) 

Here G[xi,X2, • ■ ■ , yi, 2/2, ■ ■ ■ , l/n^.; ^] is the Green's function for the hard particle fields 
with background photon field A. Ap^ and Ap^ denote the photon fields with the time 
argument being on the contours Pi and P2 respectively. 0h(p, x) is the wave function for the 
charged particle moving freely. The integration of the hard fields has two effects: 1) it gives 
the Green's function G[xi,X2, ■ ■ ■ , ?/i, 2/2, ■ ■ ■ , yuf', A]; 2) it gives radiative corrections to 
the action for the photon and electron fields, i.e. we have an effective action. Except the 
effect of renormalization, these corrections are usually small because the high energy mode 
decouples 0. 

Since G[xi, X2, ■ ■ ■ , a;„., yi, 2/2, ■ ■ ■ , yuf] A] contains only photon field, it is convenient to 
further integrate out the electron degree of freedom. Introducing source terms for the photon 
fields on paths Pi and P2, formally, we arrive at 



{R) 



n / dyfK\p'f^yf)K+^')T{ rfa:.0°(p„x.)(9^^W) 



X 



i=l 

xG 



n / rfx:0°*(p„x:)(e,W) n / dy'frM^y'f){dl+m^) 

' J=l-I 

1 5 



xG* 



xi,x2, ■ ■ ■,Xni,yi,y2, ■ ■ -^yn/, 

II III / . 

^li 2^2) ■ ■ ■ ) ^Uii yii 2/25 5 ynA ■ r T/ 



i 6 J 

where the generating functional W[J, J'] is defined by 



JW[J,J'] 



J=0,J'=0 



(3.2) 



exp{iW[J, J']} = VAV^V^ exp yiSf[A, J-A + i J' ■ A 

_ ^ji^ ^^i J dtj dxJ{t,x)-A{t,x)'^ rp ^-i J dt' J dx'J'(t',x')-A(t',x')^\^ g-j 

where we have used (• ■ ■) to denote the thermal average. Power series expansion of W[J, J'] 
can be written as a quadratic term plus higher power order terms in source fields J, J': 



W[J,J'] = -^J dxi J dx2[ J dti J rft2^^(ti,xi)J^(t2,X2)Gj,+H^i>xi;t2,X2) 
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+ 



where 



J dt[J c^4j;(t;,xOJ^(4,x2)G'(-;:H^i,xi;4,x2) 

-2 J dh J dt'Mti^^i)JUt2,^2)Q^.u{ti,^i]t'^,^2)] +■■■ (3.4) 

G^+'-\h,^ut2,^2) = (yl^(ti,xi)yl,(t2,X2))) , 

^^,(ii,Xi;i'2,X2) = (A^(ii,Xi)>l,(f^,X2)) . (3.5) 

Again, we use +, — to denote the time-ordered and anti-time-ordered product. We have 
following relations 

Gi-t'-H^i,Xi;i2,X2) = j ^e^'^-(---)^^.(t(>,<),t«,»;k) (3.6) 
where i> = max{ti, t^} and i< = min{ii, ^2} and Qf^u{t, t'; k) is defined by 

g,^(t, x; t', x') = / ^s^'^-^^'-'^'KAt, t'; k) . (3.7) 
The Green's function Qj^i, may be expressed in terms of the spectral density function: 

where PnvioJ, k) is the usual spectral density function and n{uj) is the Bose distribution factor 

- ■ (3-9) 

We shall focus on two aspects due to plasma. The first is for the momentum transfer is 
comparable with the Debye screening length, we show that the Debye screening effect can 
be properly taken into account in our formalism. The second case that we shall study is the 
modification to the transition rate due to the soft photon. 

A. Debye screening 



For the reason of simplicity, we consider the collision of two spin ^ particles and ignore 
the soft photon emission and absorption which will be considered later. We have particle 1 
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with charge qi and 2 with charge q2 moving with momentum pi and p2 initially, colliding 
with each other via EM interaction, and separating with momentum p[ and finally. When 
the momentum transfer is comparable to the Debye length for QED plasma, we shall be able 
to probe the Debye screening effect. We do a perturbation theory in QED coupling constant. 
At the leading order, 

(p'i,P2bi,P2)APi = -9i?2 / dx,e-'^^^'^'^y^^A^{x,)^l*{p[)r^l{p,) 

X / rfx2e-^(^^-P^>^^A,(x2) ^l*{p',)YHiP2) (3.10) 



where $o ^ are the wave functions of the particles 1 and 2 moving freely with specific mo- 
mentum respectively. Expression (|2.8|) for the rate reads 

ql ^Mh''<l>l{Pl) ^liPlK^M) ^{p'2h'^l{P2) ^liP2)r^l{p'2) 



(3.11) 

J=0,J'=0 



5Ux^) 6Jxix2) SJUx[) 5J',{x'2 

where VLi^2 are phase space volumes for particles 1 and 2 and the usual notation $ = (^'^'-^^ 
has been used. Since the four-point connected Green's function for photon fields is of higher 
order, we shall neglect it. Upon using Eq. ( p. 41) , 



X 



qI ^1{P\)1'^1{P1) ^liPlK^liPl) mP'2)l^^l{P2) mP2)r^l{p', 



2) 



G'^^;!{x,,X2)Gt\ ) • (3.12) 

The second term in the last square bracket corresponds to the process with particle 1 and 
2 scattering with the plasma independently. The third term is an enhancing term for the 
special case with particle 1 and 2 scattering with the plasma independently but emitting 
the photon with the same four-momenta. We shall consider only the part in the transition 
rate which involves the interaction between particle 1 and 2. This corresponds to the first 
term. Therefore, 
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><^l{P2h'^l{P2)^{p2h''^l{p'2)\{2n^^^^ (3.13) 

where k — p[ — pi is the momentum transfer and we have used 

G'^\-k) = G'^;>{k). (3.14) 

The Debye screening effect due to the plasma is taken into accounted by using the full 
two-point photon field Green's function G^~^\k). To illustrate this, we consider a simple 
case which is a collision of heavy particles, for example, protons, moving non-relativistically. 
Letting Vj-ei be the relative velocity for defining the cross section and using 

^{p[)Y^l{pi) « S,o ^l\p[)^l{pi) , ^l{p'2)r^l{p2) « S,o *0^(P'2)*0(P2) , (3.15) 
the cross section may be expressed as 

14/ ^(2-mpi+P2-p;-p'.) 

xi*j^(K)<f^(Pi)ri<f^^(P2)<f^(P2)riGi^(A;)r 

L TT^ / 7^7^(2^ S{pi+P2-Pi-P2) 



(k2+mi)) 



xi^^np;)^^(Pi)ri$^np^2)^o(P2)i %,,, , (3.16) 



where we have used to the approximated result: 

Doo(/co = 0,k)^ j^^^ (3.17) 

for k being small. Here mD is the Debye mass. When the momentum transfer is comparable 
with the inverse of the Debye screening length, it probes Debye screening length in the way 
above. 

B. Interaction between charged particle and the plasma via soft photon 



In this subsection, we shall study the effect due to the emission and absorption of the 
soft photons with energy and momentum much less than the typical momentum transfer. 
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"Soft" photons mean photons with energy and momentum much less that the momentum 
transfer. For the hard photons interacting with the charged particles, we can treat them 
perturbatively because of the smallness of the QED coupling. However, since the soft photon 
causes infrared divergences for individual diagram, we need to sum up the leading infrared 
behavior The process is always accompanied by emission and absorption of a lot of 

soft photons. Treating it perturbatively does not describe the process properly. For the soft 
photon fields, we shall sum up the leading infrared contribution to all order by using the 
Bloch-Nordsieck semi-classical method. 

The LSZ reduction formula enables us to extract the scattering matrix element from 
Green's function. Equivalently, we consider the large time behavior of the Green's function 
and find the part having the correct asymptotic wave function for the incoming and outgoing 
particles. Mathematically, 



n / dyfK*ip'f^yf)i9y^+m^)Y[ / dxi(f)l{pi,Xi){dl^+m'^)G[xi,X2, - ■ ■ ,Xn,,yi,y2, - ■ ■ ^Vuf^A] 
f=i-' 1=1'' 

= n / dy/KlPf^ yj'^ ^] n / dxi(f)h[pi, xu A]T[xi, X2, ■ ■ ■ , Xn,, yi, ?/2, ■ ■ ■ , yrif] A] (3.18) 

where (phiPi, Xi, A] and (f)h[p f , y f , A] are the wave functions of the charged particle moving in 
the classical background A field. r[a;i, X2, ■ ■ ■ , a;„^, yi, y2,---, Vnf] A] is the amputated Green's 
function for the external charged particles or the one-hard-momentum-particle-irreducible 
Green's function. Since we focus on the contribution due to the soft photons and the 
hard photon may be treated perturbatively, to the leading order, we can neglect the hard 
photon contribution. Equivalently, A carries only very small energy and momentum which 
are much less than the momentum transfer. It is then valid to neglect soft photon field 
A in r[xi, X2, ■ ■ ■ , Xn^, yi, y2,---, Vnp ^] as the leading infrared behavior is concerned since 



r describes a hard process |T0|. Our major task now is to calculate 0h[p, a;;y4] for soft 
background A field by using classical current approximation. 

For the reason of simplicity, let us first consider (p^^p^x] A\ for a scalar charged particle 
with charge q satisfying the Klein-Gordon equation: 
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(D^D'^ + m2)<^h[p,a;;A] = 0, (3.19) 

where the covariant derivative = — iqA^. It is not hard to verify that, at the leading 
order, the solution may be expressed as 

(f>^\p,x,A] = e^P^expjzg J^dz^A^'iz)'^ (3.20) 

where C is a straight line ending at with slope p^/m i.e. C is the classical trajectory of 
the charged particle moving with speed = Pfj,/m. It is generally true that the solution of 
a charged particle with any spin moving in slowly varying background A field has the form: 



Mp,x]A] = 0°(p,x)exp|zg j^dz^A^'{z)j (3.21) 
with 0[J(p, x) being the free solution since the effect due to spin is proportional to the field 



strength F^^, which is small ||TT|. It is convenient to define a classical current 

J^{z) = qv^S^^\z - X - v(^o - xo)) (3.22) 

so that we can write 

0h[p,a;; A] = 0O(p, x) expjz / dzJ^{z)A^'{z)] . (3.23) 



Now the soft photon problem can be simplified as the interaction between classical current 
and the plasma. The underlying reason for this is because the Compton wave length of the 
charged particle is much smaller than the EM field's characteristic wave length. Thus, at 
the scale of the wave length of the EM field, the charged particle behaves like a classical 
particle. This observation is originally due to Bloch and Nordsieck |l|]. 

Since T[xi,X2, ■ ■ ■ , y2, - ■ ■ , Vnp ^] is the correlation function obtained by integrat- 

ing over the hard momentum, it certainly looks local to the soft momentum. This means 
that for the purpose of treating the soft photon, we can regard it as proportional to a delta 
function for the coordinates. Let us use y = yi denote the space-time point where the hard 
collision happens. Each charged particle generates a current. The total current is the sum 
of them. Let us denote the total current as J'-^^ 
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i f 

and the Fourier transform J^j(u;,k) is simply 

i<».(.,k, ^ (e..,-^ - S^/^) . (3.25) 

where k = (cu, k). The conservation of Jj^^ is guaranteed hj J2i<li = J2f If- We now have, 
n / %'^h[P/>l//;^] n / c?a;i0h[Pi,a^i;^]r[xi,X2,---,x„^,?/i,?/2,---,l/ny;^] 

n y %<(P/> ?//) n y ^a^.^^^, x,)r[xi, ^a, ■ ■ ■ , a;„„ l/i, 2/2, ■ ■ ■ , ^n,; 0]e* i '^^■''''■^ . (3.26) 



/=1 j=l 

By translation invariance, we define F as 



h{Pi^yi)Y[ I %'/'h*(P/'2//) n / c?a;i(/){J(pi,Xi)F[xi,X2,---,x„^,?/i,2/2,---,l/nj>;0] 

/=2 j=l 

^ e^'^^^f(pl,p*2, ■ ■ ■.J^u.wLpL ■ ■ -^pO (3.27) 



where q is the total four-momentum lost for the hard particles due to the interaction with 
the plasma, i.e. 

Q = Ep'-J2p^- (3-28) 

i=i f=i 

The average rate including the correction due to the soft photon is expressed as 

(R) ^ Y,V^ n I ■ ■ ■,pW,p(,pL ■ ■ ■ / dydy'e^'^(y-y'^I{y,y') (3.29) 
where / is 



I(y,y') = e^-^'^e'''''^ (^"^^'"'X^^j,^, = e^^l-^'^''^'''''! . (3.30) 
with J(^) given by Eq. p.24| ) and J*^^'-* given by 

Jjy\z)=jjy)(z-y' + y). (3.31) 

To the leading order, we only need the quadratic term in the expansion ( |3.4D for W[J, J']. 
I{y, y') may be written as 
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/(y,t/') = exp|-^y ^ JdtJ dt'[j^{t,k)J,{t',-k)g^,{ty,t^,k) 

+j'^{t, k)Ji{t', -k)g^At^, t>, k) - 2 k)j,{t', -k)g^,{t', t, k)] ] . (3.32) 



Upon using the fact 



Jl(t,k)^J,(t-(yo-y'o),k)e^^<y-y'^ 



(3.33) 



and 



J'^{u;,k)^J,{u;,k)e^'(y'-y\ 



(3.34) 



we have 

I{y,y') ^e^pl^-tj-^^jdtj dt' [j,{t,k) - J'^{t,k)\ J,(i', -k)^^,(i', t, k)| 

(0^>^;(^'k)J.(^,k)n(^)p^.(^,k)(l-e*(^'-^))|. (3.35) 

Since current J is conserved, we can choose any gauge to work on. In Landau gauge, We 
shall use following decomposition to separate the longitudinal part and the transverse part 
of the spectral density function: 



p^,(a;, k) = p^{uj, k)P* + pi(a;, k)P; 



(3.36) 



where 



^"''^ k^- 



(3.37) 



Defining Jt,i 



jt,i = j;(cu,k)p^;;j.(cu,k) 



(3.38) 



enables to express 



Hy^y') =exp<{ - I ^ [ T^-^n{u;) {Jt pt + Ji Pi) 



2n J (27r)3 



(3.39) 
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Observing that I{y, y') depends only on y — y' in Eq. (|3.29|) , the space-time integrals over y 
and y' gives the total volume of the space-time. Therefore, 

{R) = ^ rjX^ ^l^^^i'^^, ■ ■ ■ ,pI;pLpL ■■■,piXj dye-'^ye-'^^y^ . (3.40) 
where we have defined N{y) 

N{y) ^J^J (03^(^) U Pt + Jl Pi) (l - e^'y) . (3.41) 

In principle, one should be able to use the expression ( p. 401) to compute the effect on the 
hard process due to its interaction with the plasma via soft photons. However, we shall 
not go to the detail of the computation but only focus on two aspects, namely, the infrared 
behavior and the absorption and emission of soft plasmon. Before generally discussing the 
infrared physics of formula (|3.4CI| ), we like to study first the infrared behavior of the damping 
of a charged particle moving in QED plasma with hard momentum. The charged particle 
can either move very fast or be very heavy so that its momentum is hard. 

C. Damping of a charged particle 

The question we ask is what is the probability of detecting the charged particle moving 
with the initial momentum at time Tt, i.e. the probability for finding the charged particle in 
the initial state. We expect the answer may be expressed as an exponentially decaying factor 
with the exponent proportional to the time Tt. Thus the half of the proportional factor in 
front of the time Tt in the exponent is the damping rate. We also expect that the answer may 
depend on the momentum resolution by which we judge whether the particle's momentum is 
the same as the initial momentum since the soft photon emission and absorption processes 
always appear. We shall study the leading infrared contribution to the damping of the 
charged particle due to its interaction with the plasma via the exchange of soft photons. 
Let Q, a small volume in momentum space, describe the resolution. By using the standard 
reduction technique, we can construct the probability for finding the charged particle to 
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have momentum p' lying in the phase space volume Q centered at p. Keeping the leading 
infrared terms, it may be expressed as 

P = [ ^ f rfye^(p-p')-ye-^(o.y) (3.42) 



n (27r) 



where 



= 1^1 7^M^)'^;(^'k)J.(a;,k)p,.(c.,k) (l - e^'^'^') (3.43) 



(277 



with the current 



rTt r i{uj-v-'k)Tt _ i 

J,{uj,k)= tit/rfxe'"''gV(x-vt)=gt;^— — (3.44) 

Jo J i[uj — V ■ k) 

Noting that 

J*(co', k)Jiy{uj, k) (f'v^VyTt (27r) 5{uj — v ■ k) as Tt ^ cxo (3.45) 

we can perform the integral over u and get 

iV(0, y) = q^T, j ^^(k ■ ^)v,v,p,,{k ■ v, k) (l - e^^^y) . (3.46) 

We shall look at the leading infrared behavior. The spectral density p^i, is dominated by its 
transverse part pt at the infrared region. As |k| becomes small, the leading infrared behavior 



for Pt is |T2 



Pt(^,k) 1 

— ' k2(27r)5(^) . (3.47) 

Using this, we can complete the integral over the angular part for k integral as 

JO An^k Jo ^ 



q^TTt / d(j) — -Ivl ln(l — zAIyl sin^cosi 

Jo 4:71^ 



^ |^rrt|v|ln(A|y|sin^) = |^TTt |v| ln(A|yt|) , (3.48) 

where an ultraviolet cutoff A for the k integral has been introduced, 6 is the angle between 
y and the velocity of the charged particle v, and y^ = y — (y-v)v. Therefore, 
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where Apt is the transverse size of Q with respect to the velocity v of the charged particle. 

At the leading order, the damping rate is 

g2 / A \ 

Tdamping = — T|v| In j . (3.50) 

The damping rate depends on the momentum resolution Apt logarithmically due to the 
exchange of soft magnetic photon since the magnetic force is not screened in QED plasma 
and is transverse to v. This leading piece of the damping rate vanishes if v = which is 
the case where the charged particle does not feel the magnetic force. If we take the limit 
Apt 0, the damping rate becomes infinite. This means that the particle experiences the 
long range magnetic force and can not move with a definite momentum. There has been 
discussions about the infrared problem for the damping rate of charged particles moving in 
QED plasma in the literature P]. This infrared problem can be resolved by considering only 
a physically measurable rate which introduces the momentum resolution Apt to cut off the 
infrared divergence. 

To avoid confusion, we note that only the infrared behavior of the damping rate has been 
examined. Equation ( p.50| ) is valid only when the cutoff A being larger than Ap^. Cutoff A 
is the scale below which the asymptotic behavior |3.47| can be used. At high temperature, 
A ~ eT while at low temperature A is suppressed by factor 6^^^™"= with me being the electron 
mass. Therefore, at temperature T -C rrie, this infrared cutoff becomes almost zero. Our 
resolution Apt will not be able to see this scale. Thus, the "leading infrared behavior" is 
essentially negligible. In view of Eq. ( p.48| ), this corresponds to the case A|y| ^ 1 so that 
A^(0, y) is small. These comments apply to discussions in next subsection as well. 



D. General discussion of the infrared behavior 



The cancellation of infrared divergences have been studied for the case of ideal thermal 
photon gas in reference [H. This is the case where the temperature is low and the chemical 

18 



potential for the electron is small so that the photon gas is almost ideal. Here we shall 
provide a general discussion of the infrared physics. Using the result ( |3.25| ), 



J 



t,i 



E 



E 



E 



E 



p. 



(3.51) 



Pi-k+ie J pfk—ie I \~ pi-k—ie j pj-k+ie 

Putting this into the definition ( 3.41 ) for integral N{y), it is not hard to find that there 
are singularities in the integral over u or h corresponding to pinching two poles together as 
e ^ 0. The physical meaning of these singularities which come from putting the charged 
particle's propagator on-shell is that the charged particle suffers the damping as discussed in 
the previous subsection. Therefore, their contributions to N{y) can be expressed similarly as 
the result shown in Eq. (|3.48| ). We make following replacement to subtract out the damping 
effect of the single particle: 



j;(k)Uk)-.j;ik)Mk)-Y: 



{pi - k + ie){pi ■ k-ie) 



E 



{pf ■ k - ie){pf ■ k + 



(3.52) 



Because of Debye screening, the exchange of the magnetic photons dominate at the infrared 



region for the u, k integrals.. By the sum rule |[T2 

du Pt{uj,k) 



27r UJ 



1 

k2 



(3.53) 



and the behavior 



Pt('^,k) 



UJ 



j^(27r)5(u;) 



(3.54) 



as k — s> 0, we have the leading infrared contribution to N{y) 



infra 



T 



d^k 1 
(27r)3k2 



J*(0,k)-J(0,k)-^ 



Vj ■ k+ie 



E 



Vf ■ k—ie 



(3.55) 



besides the contribution from the damping of the individual particles. The k integral in 
expression above may be written as a one-parameter integral as shown in the appendix. 
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Using the conservation of the current so as to replace P*^ by 5ij , apart from the contribution 
due to the damping effect, 



Niy) 



infra 



T\y\ 



i=l i'=i+l f=l /'=/+! 



ni rif 



i=i f=i 



(3.56) 



where J(y,v,v') is defined and simphfied in the appendix. Since J(y, v,v') contains no 
infrared divergence, A^(y)mfra is finite. The unscreened static magnetic force in QED plasma 
causes no infrared divergences. 



E. long wave plasma oscillation 

Besides to the leading infrared effect due to the unscreened static magnetic force, there 
are also effects due to the excitations of the long wave plasma oscillation. To simplify the 
discussion, let us isolate the contribution of the plasmon by considering 

pt,i = (27r)5(a;2-u;t,i(k)2) (3.57) 

with u = cut.Kk) being the dispersion relation for the transverse and longitudinal plamons 
respectively. Putting this into the expression (|3.41|) gives 



+ 



Ji 



+ 2n(u;t(k))(l - coskty) 



"(1 - e*^'^) + 2n(cui(k))(l - cos%) 



(3.58) 



2^1 (k) 

where fct,i = (i^t,i(k), k). For the integrals in Eq. ( |3.58| ), the plasma mass mpi = c<;t,i(k = 0) 
cuts off the integral at the infrared region. At very high temperature rripi ~ eT or at 
low temperature with sufficiently large chemical potential, rripi ^ e^T, N{y) is of order 
e^T/rripi <^ 1, thus perturbation theory works well since we do not need to sum over the 
series to form the small exponential factor. 

At low temperature and small chemical potential, rupi may be much smaller than e^T. 
There are two cases that we would like to discuss. 
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1) For the case where the exchange of energy between the charged particles and the 
plasma is much smaller or the same order of magnitude of rripi, we then need to compute 
integral ( |3.58| ) by using, for example, the one-loop spectral density function. The calculation 
is rather involved. We leave it to another day. 

2) For the case where the exchange of energy between the charged particles and the 
plasma is much larger than mpi, we can set rripi to be zero and simply use the dispersion 
relation co't,i(k) = |k|. For this case, the contribution to N{y) from the plasmon dominates 
since the photons are nearly ideal. Therefore, 

d^k 1 



N{y) ^ iV(y)pi 
Writing 



(27r)3 2|k| 



M\k\,k)\' (1 - e^'y) + 2n(|k|)(l - cos%) 



(3.59) 



l^.(|k|,k)|^ = 
with four- vector m = (1, k), we have 



1 

k2 



rii 



(3.60) 



N{y) 



dk 



QfPff^ 



j^i Pf^ 



doj 
2uj 



(1 



+ 



[1— cos Lou-y) 



(3.61) 



Expression (|3.61| ) above contains an ultraviolet divergence which should be cut off at the 
hard momentum scale A. Instead of using a sharp cut off, we use an exponential factor 
g-^w/A suppress the contribution from the region with u bigger than A. Using the result 
in the appendix for the u integral. 



N{y) 



1 



167r3 



dk 



hi Pi-u f^i PfU 



In (l+iAu-y) + In 



sinh TrTu-y 
TcTu-y 



(3.62) 



A full analysis of expression ( p.62|) is still formidable. We shall consider a special case 
where the process involves at least one heavy particle moving non-relativistically so that 
it can absorb the extra momentum without affect much the conservation of energy. Thus, 
we effectively has only the energy conservation without worrying about where to dump the 
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extra momentum. For this case, we can effectively replace y by and only left with the 
time component. It is convenient to consider the differential rate with the energy exchange 
with the plasmons in plasma to be in the interval {E, E + dE): 

sinh TiTt 



d{R) oc J dte-'^^expl - A 



ln(l + iAt) +ln' 



TlTt 



dE 
'2^ 



dE 
2^ 



(3.63) 



where A is defined by 



A 



9/ 



(3.64) 



IGn-^ J \ " " Pi-u t^-, " pfU 
We can then reproduce two previous results For the case T ^ i?, it is expected 

that the result is approximately the same as the zero temperature case. For the t integral 
in expression ( p.63| ), typical t is of scale 1/E and thus Tt <^1. Noting sinh a; ^ x for small 
X, we have 



/d F 
dte-'^\l + iAt)-^ — ~ dE9{E) 
2tt 



dE A 
27r ''~'E \a) 



(3.65) 



Therefore, the transition rate with soft photon carrying away energy no more than AE with 
> T is 



f 

{R)ae oc / 
Jo 



A /AE\^ 



dE - [- 
E VA 



(3.66) 



This is the zero temperature result 0. 

For the case where the energy exchange between plasma and charged particles are much 
smaller than the temperature, T ^ E. Now for typical t in the integral in Eq. ( p. 63 ), 
Tt^ 1. Using sinhx ~ for large x. 



dE /T\ 

27 va; 



AJ 



dE 



TiAT+iE TxAT-iE 
2'kAT 



27r ^2 + (vrAT)2 ' 



(3.67) 



The agrees with previous result 
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IV. CONCLUSIONS 



A general formula for computing the transition rate for hard process happening in a 
soft thermal background is derived. In particular, we use it to study hard process involving 
charged particle inside QED plasma. Both the Debye screening and the damping of the 
hard charged particles are incorporated in our formula. We show that the transition rate is 
free of infrared divergence even though the magnetic force is not screened for QED plasma. 
Soft plasmon emission and absorption are analyzed. Our formalism can also be applied to 
the case where some or all of the hard momentum particles are replaced by particles which 
interact weakly with the plasma so that after the formation of the plasma, it does not have 
enough time to have thermal contact with the degree of freedom corresponding to these 
particles. 
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APPENDIX A: COMPUTATION OF INTEGRALS 



We first evaluate following integral: 

iy,v,v;_ J (27r)3vk + iev'-k + iek2 



Exponentiating the denominators produces 

I{y,v,V)^--- dt dt' I - 

\y\ Jo Jo 

Completing the k integral gives 



(27r)3 k2 



,ik-(vt+v't') 



2 - e^^-y - e 



ik-y 



1 roo roo 

/(y,v,vO = -^ / dt/ dt' 
y -^0 Jo 



\vt + v't'\ \vt + v't' + y\ \vt + v't'-y\ 



(Al) 



(A2) 



(A3) 



We can insert another integral over 5 and change the variables t — sa^ t' — sf3 to get 



poo fj Q roo roo 

7(y,v,v') = -/ j^8{s-t-t') dt dt' 
Jo \y\ Jo Jo 

^0 y ^0 



\vt+v't'\ \vt+v'f+y\ \vt+v't'-y\ 
1 1 



av+{l—a)v'\ \av+{l—a)Y'+y/s\ |q;v+(1— q;)v'— y/s| 



lo '^'^[av+(l-tt)v']^ lo 
Jo 



ds 



2 + cos 9{a) In 



^l-2scos^(a)+s2 ^1+2scos^(q;)+s2 
1 + cos 6'(q;) 



1 — cos^(q;) 



[av+(l-«)vf 

where 9{a) is the angle between vectors q;v+(1— q;)v' and y. 
We now evaluate integral 

~ 2 



(A4) 



Jo 

Expanding the denominator gives 



X — 1 



[1 — COSTjx) . 



n=l 



dx 

X 



2e 



g nx+ir)x ^—nx—ir]x 



1 + 



In 



n=l 



sinh 717] 
Tirj 



where we have used 



r°° dx 
Jo X 



e-^^) = In ^ 



and the infinite product 



sinhx = 3; I H- 



X 



n=l 



(A5) 



(A6) 



(A7) 



(A8) 
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